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Summary 

This report investigates the feasibility of using a digital filter to change the pre- 
emphasis applied to an audio signal encoded in digital form from the CCIR (50 p.s) 
characteristic to the CCITT characteristic. This filter is required at the interface between 
bit-streams carrying 13-bit and companded 10-bit digital sound signals. The disadvantage 
of using an analogue method for this purpose is that a separate filter would be needed for 
each audio channel in the multiplexed signal, preceded by a digital-to-analogue converter 
and followed by an analogue- to-digital converter. 

It is shown that satisfactory results would be obtained by using a relatively 
simple (second-order) digital filter, and with filter coefficients specified to 8-bit accuracy. 
Some hardware aspects are considered in the report, and it is concluded that, by using 
suitable time sharing techniques and appropriate data storage, one central (multiplier) 
unit could serve a group of six multiplexed audio channels, which would be convenient 
in one application of the companded digital sound system. 
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DIGITAL SOUND: FEASIBILITY OF CHANGING FREQUENCY 

PRE-EMPHASIS BY A DIGITAL FILTER 

T.A. Moore, B.E., M.Eng.Sc. 



1. Introduction 

A sampled-data filter consists of an interconnection of 
delays, multipliers and adders designed to implement a 
specific frequency shaping operation on a signal; a digital 
filter is a sampled-data filter in which the input signal has 
been digitised and in which the processing is carried out 
digitally. 

The work described in this report was undertaken to 
assess the feasibility of using a digital filter to alter the pre- 
emphasis applied to audio signals from the CCIR (50,us) to 
the CCITT pre-emphasis characteristics. 1 This filter is 
intended for use at the interface between a 6-336 Mbs" 1 
bitstream carrying thirteen 13-bit linear sound channels and 
a 2-048 Mbs -1 bit stream carrying six 10-bit companded 
sound channels. 2 However, for the purposes of this report 
the companded system may be regarded as a 13-bit linear 
system using a different pre-emphasis, because the system 
uses a 13-bit linear signal as a starting point for digital 
companding in such a way that the quantising noise is sub- 
jectively close to that of a 13- bit linear system. The 
sampling frequencies for these sound signals would be 
32 kHz. The alternative to using a digital filter would be 
to employ an analogue frequency shaping network pre- 
ceded by a digital-to-analogue converter (d.a.c.) and 
followed by an analogue- to-digital converter (a.d.c). Apart 
from the usual advantages of digital operation, there is the 
possibility of 'sharing' one digital filter between a number 
of sound channels in the multiplexed signal, thereby re- 
placing by one unit that number of analogue conversion 
units. 



The desired transfer characteristic is shown in Fig. 1 : 
it is given in explicit form by the equation: 
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where H(jco) is the transfer function corresponding to an 
angular frequency co = 2-nf (where /is the frequency in Hz). 
The response is specified from 40 Hz to 15 kHz. (The gain 
factor, 2-279 is chosen to give a maximum gain of unity.) 

This frequency characteristic may be realised by a 
simple (second-order)* analogue network. In certain 
instances a simple form of design, based on modifications 
to the corresponding analogue filter design theory, yields a 
digital filter of about the same complexity as that of the 
equivalent analogue filter; 3,4 however this requires that 
either the amplitude-frequency characteristic be constant 
within given frequency segments (as would be the case with 
a low-pass, high-pass or band-pass filter) or that the highest 
frequency of interest be low compared to half-sampling 
frequency. The present requirement does not fall into 
either of these categories, and therefore a more complicated 
technique is needed to design the digital filter, as shown in 
Section 3. 

* The transfer function, G(jco), of an analogue filter can generally 
be written as a ratio of two polynomials in jco; the order of the 
filter is then the highest power of joj in the denominator of G(jCO>. 
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Fig. 1 - 50 jds to CCITT pre-emphasis transfer characteristic 
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output 



Fig, 2- General form of linear recursive sampted-data filter 
I a n = A Q a n torn = 1, 2 . . . m 

indicates one clock-period delay (ail elements clocked simultaneously) 



2. Theory 

This section outlines the form of a digital filter, and its 
behaviour. 

If the input signal is sampled, then a device which 
operates in some fashion on those sample values is called 
a sampled-data filter. If the sample values are represented 
digitally and the operations (such as addition and multi- 
plication) within the filter are carried out using digital 
methods, then the filter is digital. 

2.1. Sampled-data filters 

In a sampled-data filter the input signal is sampled {at 
a regular rate, 1/T) and the filter output, produced at each 
sampling instant is a weighted sum of the input, several past 
(i.e. delayed) inputs and several past outputs. If the output 
does not depend on previous values of the output, the filter 
is non-recursive and has a finite-length impulse response; 
otherwise the filter is recursive and has an infinite impulse 
response. 

The number of past input or output values (whichever 
is the greater) determining the filter output is termed the 
order of the filter; obviously the filter complexity increases 
with its order. 

The form of a sampled-data filter is shown in Fig. 2. 
The equation describing the behaviour of the filter Is called 
the difference-equation of the filter and is given by 



-b x y k _ x -b 2 y k _ 2 - 



(2) 



where u k and y^ are the filter input and output respectively 
at the fcth sampling instant. 

A frequency response of the filter is associated with 
each set of coefficients (A Q> a v a v . . . . b { , b 2 ....). 
Denote the Fourier Transform of the sequence u k by U(jgj) 
and ofy^ by Y(jaj). Then by taking the Fourier transform 



of Equation (2), and noting that if a time function is delayed 
by time T its transform becomes multiplied by e~ >ojT , the 
frequency response G(jto) of the filter is given by:— 

Y(jw) 1+a ie ~ ja;T +^e- |2wT + ... 
G(M=— ~- = A n . ,' _ ;/ „ T 2 _ ; „, <T (3) 



U(jw) 



■\+b l e-' i " T +b 2 e-i 20JJ + ... 



G(jco) is a complex number, whose modulus and 
phase at a given value of o> yield the amplitude and phase 
responses of the filter at that (angular) frequency. 

By its nature the frequency response is periodic with 
frequency, the period being the sampling frequency. How- 
ever, the portion of the response relevant to most appli- 
cations is that in the region below half-sampling fre- 
quency. 

2.2. Digital filters 

Because the coefficient values as well as the signal 
values in a digital filter are represented digitally, in any 
practical situation quantisation will occur and the filter 

behaviour will differ from that of a corresponding sampled- 
data filter. 

The length of the digital word representing each 
coefficient has a bearing on the complexity of the filter, and 
in any design is kept as short as possible while maintaining 
the filter response within the desired specification. The 
extent to which each coefficient may be allowed to deviate 
from its ideal (i.e. analogue) value depends on the sensitivity 
of the filter characteristic to parameter changes at the 
critical points on the response. 

Also inherent in a digital filter is data (or signal) 
quantisation. In the present application, the input and 
output signals are specified to 13-bit accuracy. Certain 
products within the filter will have a higher accuracy: for 
example if N-bit coefficients are used, some products com- 
prise (13 + N) bits. There must, however, be a truncation 
at the filter output if only the 13 most significant bits 
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describing the signal are used for transmission;* this is 
equivalent to re-quantising the signal and means that the 
quantising noise power at the digital filter is (roughly) 
twice that at the input In this respect a digital filter 
would not be superior tc an analogue filter with associated 
d.a.c.'s and a.d.c.'s. 

In order to minimise errors due to truncation, the 
values of .yjt-v y~k-2 snou 'd be specified to more than 13- 
bit accuracy before they are multiplied by b v b 2 : in other 
words the truncation to 13 bits at the filter output is a 
'final' truncation and is not applied to those values fed back 
to the filter, such as^-v ^jfc-2- However, it is necessary 
to apply some truncation in these paths since if this were 
not done the number of bits in the untruncated outputs 
would increase by N at each sampling instant. 

In a well designed filter the quantising noise power 
introduced by this internal truncation is well below the in- 
put (and output) quantising noise power. The effect of 
this noise on the filter output signal, in magnitude and 
spectral content, depends on the particular configuration in 
which the filter is realised; the amount of noise produced 
at the output may be estimated by imagining a white noise 
source ateach internal truncation and by tracing its progress 
through the filter. (See Appendix.) 



3. Filter design 

Because of the difference in the forms of the frequency 
responses of a digital and analogue filter, any digital filter 
designed to simulate the response of an analogue filter 
would only do so approximately; this approximation 
would improve as the order of the digital filter was 
increased. 

An estimate of the allowable deviation between the 
two responses for a filter intended to change the pre- 
emphasis of a sound signal from the CCIR to the CCITT 
characteristic may be judged from Table 1, the specification 
of the frequency response tolerance of an audio link. 

TABLE 1 

Frequency Response Tolerance of an Audio Link 



Frequency range 

40 Hz to 125 Hz 
125 Hz to 10 kHz 
10 kHz to 14 kHz 
14 kHz to 15 kHz 



Tolerance 

+0-5dBto-2dB 

+0-5 dB to-0-5dB 
+0-5dBto-2dB 
+0-5 dB to -3 dB 



This section outlines the design of the digital filter for 
the type of application mentioned above. Most of the 
effort concentrated on trying to find whether a simple (i.e. 
low-order) sampled-data filter could be designed with an 



Although the present application is a 10-bit companded signal, the 
latter signal may be taken as a 13-bit linearly coded signal as 
explained in Section 1. t 



acceptable approximation error. Following this, the effect 
of quantising the filter coefficients was investigated, as were 
the data-truncation effects within the filter. Finally, some 
thought was given to the hardware aspects of the filter. 

3.1. Design of sampled-data filter 

3.1.1. Invariant impulse response method 

When designing a sampled-data filter to simulate the 
response of an arbitrary analogue filter, a standard technique 
if that of invariant impulse response. As its name suggests, 
this technique yields a sampled-data filter whose impulse 
response matches the sampled impulse response of the 
analogue filter; the order of the digital filter is then equal 
to that of the analogue filter. However, this technique can 
only be applied if the response of the analogue filter is 
band-limited to half-sampling frequency. 

In the present application the (second-order) analogue 
filter has no inherent low-pass characteristics and therefore 
a (conceptual) low-pass filter is required in tandem to 
attenuate the continuation of the filter response at fre- 
quencies above half-sampling. Assuming a 32 kHz sampling 
frequency, this low-pass filter should have a passband from 
to 15 kHz, a transition band between the pass-band and 
stop-band 2 kHz wide, and also the allowable deviation in 
the pass-band should not exceed about ±0-1 dB. The order 
of such a low-pass filter is estimated as 4. This gives the 
combined filter an order of 6, which is high; the corres- 
ponding sampled-data filter would thus have an order of 6. 
In addition, employing this technique for the present 
application would yield an over-designed filter, in the sense 
that its phase response would be matched to that of the 
analogue filter; this is unnecessary for sound signals. This 
design method was therefore discarded in favour of the 
method described below. 

3.1.2. Empirical method of design 

An empirical approach was next taken in designing 
a suitable sampled-data filter. In this case the parameters 
(i.e. the coefficient values) of a recursive digital filter of a 
given complexity (order) are optimised, generally by means 
of an algorithm, to obtain a best fit between the analogue 
filter response and the sampled-data response given by 
Equation (3). 

An optimisation procedure was carried out on the 
computer using the Fletcher-Powell algorithm, the error 
function being minimised was the sum of the squared dif- 
ferences, in decibels, between the analogue and digital filter 
amplitude responses at a number of frequencies uniformly 
spaced over the frequency band of interest (40 Hz to 15 
kHz). The algorithm was run using 30 points on the fre- 
quency scale; it yielded a set of coefficient values for a 
second-order filter which produced an approximation error 
of less than 0-25 dB across the band. 

3.1.3. Adjustment of parameter values 

The next stage in the design was to vary the para- 
meter values of this second-order sampled-data filter slightly 
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to see whether the maximum approximation error could be 
reduced in those regions where the tolerance was most 
stringent (i.e. 125 Hz to 10 kHz) at the expense of a greater 
error outside that region. By evaluating the sensitivity of 
the {log} frequency response to changes in the filter para- 
meters it was possible to judge the effect of varying the 
parameters slightly. In the case of the second-order 
sampled-data filter the following set of coefficients (using 
the notation of Equations (2) and (3)} was arrived at. 



A. 



a. 



a„ 



'0 "1 ~2 1 2 

0-661 -0-6695 -0-2286 -0-9050 -0-1174 



be coded to 13-bit accuracy. The effect of reducing the 
coefficient word lengths was next investigated, this being 
done by setting each coefficient value to the nearest integer 
multiple of 1/1024 and evaluating the resultant filter fre- 
quency response. By again exploring the sensitivity of the 
(log) frequency response to changes in parameter values {as 
in Section 3.1) it was possible to estimate the effect of 
further quantising the coefficients; these results were used 
to make the multiples of 1/1024 divisible by as high a 
power of 2 as possible, while maintaining the approximation 
error within acceptable limits. The following set was 
arrived at: 



The corresponding approximation error versus fre- 
quency* is shown by curve {a) in Fig. 2. It can be seen 
that the error at any frequency is less than 20% of the 
'allowable' tolerance for an audio link quoted in Table 1. 

Since the approximation error was judged to be 

acceptable, the parameter values quoted above were taken 
as the basis of a second order digital filter designed to 
match the analogue response of Equation (1). 

3.2. Quantising effects 

3.2.1 . Coefficient word length 

To represent each coefficient digitally to the 
accuracy of the above values, each coefficient would need to 



The sense is such that the error is positive where the digital fitter 
gain exceeds the desired gain. 



677 



-684 



-228 



-928 



188 



1024 1024 1024 1024 1024 

(0-6611) (-0-6680) (-0-2227) (-0-9062) (0-1836) 

The corresponding approximation error is shown by curve 
{b) in Fig. 3; this curve deviates from the ideal response by 
less than a quarter of the tolerances quoted in Table 1. For 
each coefficient the numerator is divisible by 4 (except for 
the gain parameter A Q ), and therefore they may all be 
specified by 8 bits. 

3.2.2. Filter structure 

The effect of signal quantising noise (or data- 
truncation effects) depends on the digital filter 'structure', 
or configuration. Two such structures are shown in Fig. 4; 
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Fig. 3 - Approximation error vs. frequency 

(u) Unquantised coefficients ib) Coefficients quantised to 1 part in 256 
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{a) 




(b) 



Fig. 4 - Second-order digital filter structures 
indicates one clock-period delay (all elements clocked 
simultaneously) 



these are considered because, as required for experimental 
reasons, it is desired to produce a filter whose gain is 
dependent on only one parameter, and to keep any such 
variable parameter on the output side of the filter. Com- 
paring those structures with that of Fig. 2 it will be noted 
that the number of word stores is halved by suitable arrange- 
ment of the filter constituents. 

In Fig. 4, N x and N 2 represent the quantising noise 
introduced due to internal truncation of the products 
b } y^_i and b^y^_2. The Appendix outlines how the 
spectrum of the quantising noise in the output signal due 
toiVj and-/V" 3 may be calculated. 

For the parameter values mentioned above, there is 
not much to choose between {a) and {b); however (b) is 
preferred because the same spectrum shaping would be 
applied to the truncation 'noise' as to the quantising noise 
embedded within the input signal {«). 

3.2.3. Signal word lengths 

A good choice for the store word-length (of Fig. 
4{a)) would be 20 bits. For a 13-bit input signal, this 
would allow the most significant bits (m.s.b.) of the stored 
words to be 2-bits (or 4 times) more significant than the 
m.s.b. of the (13-bit) input. The 2 extra significant bits 
would allow for the fact that the stored word is obtained 
by adding together 3 separate numbers; the 5 least signifi- 
cant bits of the stored words would render the r.m.s. quan- 
tising noises N { and N 2 to be 2 -s times the input signal 



quantising noise. The effect ofiVj and N would therefore 
be to increase the overall signal-to-quantising noise by at 
most 0-005 dB; this would be a negligible addition to the 
3dB increase in quantising noise due to the 'final' truncation 
to 13 bits discussed earlier, 

3.3. Hardware considerations 

From the discussion in the previous paragraph, 20 
bits appears to be over-generous as a specification for the 
stored word lengths; however it is convenient to have word 
lengths as multiples of 4, and a word length of 16 would not 
be quite enough. 

The most complex item in the filter is the multiplier 
unit, used to multiply different signal combinations by the 
coefficients A Q , a v a v b x and b 2 . With currently available 
devices it is possible to multiply a 20-bit by an 8- bit 
number and add the product to a running sum in 700 ns. 
By suitable storage and timing arrangements the five multi- 
plications necessary for a second-order digital filter could 
be carried out in a single multiplier on a time-shared basis in 
about 4 ^s. Since the sample period for digital sound is 
31-25jUs, the filter could be time-shared between 6channels. 
This form of filter would require the storage of 12 words 
(i.e. 2 for each channel comprising the multiplexed audio 
signal), each 20 bits long. The amount of shift-register 
hardware might be reduced by using serial storage instead 
of parallel (i.e. for the parameters quoted, using two 120- 
bit instead of 40 6- bit registers); however, this would 
necessitate the up-conversion of the clock rate from 32 kHz 
at certain points within the filter. 

On current prices, the total cost of the digital filter 
components is estimated to lie in the region £100 to £200. 
This compares favourably with the estimated component 
costs for analogue filtering (including a d.a.c. and a.d.c. for 
each of the six channels) of about £1,500; the costs 
exclude those of the associated demultiplexer and output 
multiplexer required in each case. 



4. Conclusions and recommendations 

It has been shown that the digital conversion of audio 
signals from the CCIR (50 us) pre-emphasis characteristic 
to the CCITT frequency characteristic is feasible using a 
simple (second-order) filter whose coefficients are specified 
to 8-bit accuracy. Because of the possibility of time- 
sharing the filter between six audio channels and of using 
one multiplier 'unit' on a time-shared basis to carry out all 
the multiplier operations required for a filtering operation 
the amount of hardware in the filter would be quite modest. 

It is important to note the overall degradation in the 
signal-to-quantising noise ratio is similar whether digital or 
analogue coupling is used. However, on basic component 
costs, the digital filter would have a significant advantage 
over a set of d.a.c. — analogue filter — a.d.c. combinations, 
one of which would be required for each audio channel if 
analogue coupling were employed. Therefore, it is recom- 
mended that serious consideration be given to the develop- 
ment of a digital conversion system for use at the interface 
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between the 13-bit linear and 10-bit companded digital 
systems, the deciding factor being the extent to which such 
a conversion would be likely to be used during a period 
when both systems would be in operation. The choice 
would be largely economic, but the added reliability of the 
digital method should also be taken into account. 
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6. Appendix 
Analysis of Internal Data Truncation Effects 



It is convenient to employ z-transform notation to 
calculate the effect of the quantising noises,// of Figs. A{a) 
and 4(b), on the output signal. 

The z-transform of a sequence juA is defined as 




u k-]*] 




U(z) = 



the Fourier Transform of a series of impulses amplitude 
modulated by this sequence 




w k 6(t-kT) 



k=-< 



i.e. the output is the discrete convolution of the input 
sequence and filter impulse response. 

Since the z-transform of a sequence x k = «k_i is 
X(z) = z _1 U(z), the z-Transfer Function of single delay is 



Referring to the digital filter structures of Fig. 4, the 
quantising noises (N), due to data truncation effects are 
assumed to be independent white noise sources, of zero 
mean value and of power a 2 /12 where a is the quantising 
step of the truncated products. 

The equations governing the filter of Fig. 4{a) are: 



is then 



U(e 



jOJT. 




UuQ 



-jkCOT 



k=-c 



V'(z) = z" 1 [a 2 U(z) + N;<z) - b 2 X'(z)] 

W'(z) = z" 1 [a t U(z) + N'j (z) - Z>,X'{z) + V'(z)l 

X'{z) = W'(z) + U'(z) 



Y'{z)=A X'(z) 



T being the sample period. 

Using this notation the Transfer Function G(z) of a 
sampled-data filter is the z-transform of its Impulse 
Response IgA which is the response to an input u Q = 1, 
u k = for k ¥= 0; this means that if the system has any 
input sequence \u k l the z-transform of the output sequence 
\y k ] is obtained by multiplying the system transfer function 
by the z-transform of the input sequence; thus 

Y{z)=G(z}U(z) 

In the 'time' domain, 



Re-arranging: 

Y'(z) = G(z) U(z) + G'j (z)N', (z) + G; (z)N; (z) 



where 



G{z)=^l 



1 +«,z" 



+ £l 2 2 



g;<z) 



l+^z" 1 +b 2 z~ 2 



z = 



V" 1 



^+b l z~~' +2> 2 z -2 



V" 2 



1 +6,z _l +b 2 z~ 2 
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The portions of Y'(z) due to N\ and N' 2 thus have 
spectral densities of 



a' 
12 



G>I WT ) 



2 a 
and— Gi{e ,coT J 
12 2 



respectively. 



g;'(z) = G(z) 



The spectral density of the truncation noise at the 
output of the filter is then 



The equations governing the filter of Fig, 4(6) are: 
V"(z)=U{z) + N' 1 '(z)+N;'( Z )-& 1 z- 1 V"{z)-i 2 z- 2 V"{ Z ) 
Y"(z)=^ [V"(z) + fl 1 z- 1 V"{z) + a 2 z~ 2 V"(z)] 

Re-arranging: 

Y"(z) = G(z)U(z) + G';(z)N';(z) + G' 2 '(z)N;'(z) 
where G"(z)=G(z) 



a- 
12 



G(e 



JWTi 



, for N j and N, 



Note that due to the form of the filters in Fig. 4, it is 
not strictly correct to say that N t and N 2 arise due to the 
truncation of b 1 y k __^, b 2 y k _2', the products being 
truncated are actually 



b x y 



k-1 



Vk-2 



and 



A. 
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